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Solar radiation is a primary driver for many physical, chemical, and biological processes on the earth's 
surface. Solar energy engineers, architects, agriculturists, hydrologists, etc. often require a reasonably 
accurate knowledge of the availability of the solar resource for their relevant applications at their local. 
In solar applications, one of the most important parameters needed is the long-term average daily 
global irradiation. For regions where no actual measured values are available, a common practice is to 
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estimate average daily global solar radiation using appropriate empirical correlations based on the 
measured relevant data at those locations. These correlations estimate the values of global solar 
radiation for a region of interest from more readily available meteorological, climatological, and 
geographical parameters. The main objective of this study is to chronologically collect and review the 
extensive global solar radiation models available in the literature and to classify them into four 
categories, i.e., sunshine-based, cloud-based, temperature-based, and other meteorological parameter- 


TM based models, based on the employed meteorological parameters as model input. 

Furthermore, in order to evaluate the accuracy and applicability of the models reported in this paper 
for computing the monthly average daily global solar radiation on a horizontal surface, the geographical 
and meteorological data of Yazd city, Iran was used. The developed models were then evaluated and 
compared on the basis of statistical error indices and the most accurate model was chosen in each 
category. Results revealed that all the proposed correlations have a good estimation of the monthly 
average daily global solar radiation on a horizontal surface in Yazd city, however, the El-Metwally 
sunshine-based model predicts the monthly averaged global solar radiation with a higher accuracy. 
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1. Introduction 


Solar radiation arriving on earth is the most fundamental 
renewable energy source in nature. Solar energy, radiant light 
and heat from the sun, has been harnessed by humans since 
ancient times using a range of ever evolving technologies. The 
energy from the sun could play a key role in de-carbonizing the 
global economy alongside improvements in energy efficiency and 
imposing costs on greenhouse gas emitters. In the studies of solar 
energy, data on solar radiation and its components at a given 
location are a fundamental input for solar energy applications 
such as photovoltaic, solar-thermal systems, solar furnaces and 
passive solar design. The data should be reliable and readily 
available for design, optimization and performance evaluation of 
solar systems for any particular location. The best way to 
determine the amount of global solar radiation at any site is to 
install measuring instruments such as pyranometer and pyrheli- 
ometer at that particular place and to monitor and record its 
day-to-day recording, which is really a very tedious and costly 
exercise [1]. In spite of the importance of solar radiation mea- 
surements, in many developing countries this information is not 
readily available because of not being able to afford the measur- 
ing equipments and techniques involved. Therefore, a number of 
correlations and methods have been developed to estimate daily 
or monthly global solar radiation on the basis of the more readily 
available meteorological data at a majority of weather stations. 
Empirical models which have been used to calculate solar radia- 
tion are usually based on the following factors [2]: 


1) Astronomical factors (solar constant, earth-sun distance, solar 
declination and hour angle). 

2) Geographical factors (latitude, longitude and elevation of the 
site). 

3) Geometrical factors (azimuth angle of the surface, tilt angle of 
the surface, sun elevation angle, sun azimuth angle). 

4) Physical factors (scattering of air molecules, water vapor content, 
scattering of dust and other atmospheric constituents such as O>, 
N2, CO», O, etc.). 

5) Meteorological factors (extraterrestrial solar radiation, sun- 
shine duration, temperature, precipitation, relative humidity, 
effects of cloudiness, soil temperature, evaporation, reflection 
of the environs, etc.). 


The number of correlation that have been published and tested 
to estimate the monthly average daily global solar radiation is 
relatively high, which makes it difficult to choose the most 
appropriate method for a particular purpose and site. Selecting 
an appropriate methods from various existing models is based on 
their data requirements (the selected methods utilize only daily 
variables normally available at a majority of weather stations) 
and the model accuracy. 

Estimations of the monthly mean daily global solar radiation 
for a large number of locations are presented in various research 
works. The main objective of this study is to comprehensively 
collect and present the global solar radiation models available in 
the literature, including the study carried out on the estimation of 
the monthly average daily global solar radiation on horizontal 
surfaces and to classify them into four categories based on the 


employed meteorological parameters as model inputs. It is 
noteworthy that all the proposed models contain empirical 
constants which depend on the season and the geographical 
location of a particular place. The empirical coefficients for some 
correlations available in the literature are presented. The models 
collected and reviewed in this research are chronologically pre- 
sented and therefore are useful for selecting the appropriate 
model to estimate global solar radiation for a particular place of 
interest. 

In addition, in order to evaluate the performance of the models 
in each of the four categories described in Section 2, some 
selected models from each category are used to estimate the 
monthly average daily global solar radiation by using, astronom- 
ical, geographical, geometrical and meteorological data of Yazd 
city, Iran. The solar radiation values produced by the selected 
models are then compared with each other as well as with the 
available experimental data. Among the selected models in each 
category, the model which produces the most accurate values of 
global solar radiation is selected and used for predicting the 
monthly average daily global solar radiation in Yazd. 


2. Modeling of global solar radiation 


Solar researchers have developed many empirical correlations 
which determine the relation between solar radiation and various 
meteorological parameters. As the availability of meteorological 
parameters, which are used as the input of radiation models is the 
most important key to choose the proper radiation models at any 
location, empirical models can be mainly classified into four following 
categories based on the employed meteorological parameters: 


(1) Sunshine-based models. 

(2) Cloud-based models. 

(3) Temperature-based models. 

(4) Other meteorological parameter-based models. 


Among all such meteorological parameters, bright sunshine 
hours, cloud cover and temperature are the most widely and 
commonly used ones to predict global solar radiation and its 
components at any location of interest. 


2.1. Sunshine-based models 


The most commonly used parameter for estimating global 
solar radiation is sunshine duration. Sunshine duration can be 
easily and reliably measured, and data are widely available at the 
weather stations. Most of the models for estimating solar radia- 
tion that appear in the literature only use sunshine ratio (S/So) for 
prediction of monthly average daily global radiation. In the 
following section, correlations which use only the sunshine ratio 
as the key input parameter are presented and classified based on 
their developing year. 


2.1.1. Model 1: Angstróm-Prescott model 

The first and the most widely used correlation for estimating 
monthly average daily global solar radiation was proposed by 
Angstróm [3], who derived a linear relationship between the ratio 
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Nomenclature 

a-t regression coefficients 

C mean total cloud cover during daytime observations 
(octa) 

CF cloud factor 

E mean evaporation (cm) 

H monthly average daily global radiation on horizontal 
surface (MJ/m? day) 

Ho monthly average daily extraterrestrial radiation on 
horizontal surface (MJ/m? day) 

Hc monthly average clear sky daily global radiation on 
horizontal surface (MJ/m? day) 

Hp monthly average daily beam radiation on horizontal 
surface (MJ/m? day) 

Ha monthly average daily diffuse radiation on horizontal 


surface (MJ/m? day) 

the estimated global solar radiation corrected for 
systematic bias (MJ/m? day) 

Haay half-day length 


I, hourly beam radiation (MJ/m? h) 

Iq hourly diffuse radiation (MJ/m? h) 

Isc solar constant (1367 W/m?) 

n number of day of year starting from the first of 
January 

P mean precipitation 

Po atmospheric pressure at sea level (hpa) 

Ds atmospheric pressure at the site (hpa) 

PS ratio between mean sea level pressure and mean 


daily vapor pressure 


PWV precipitable water vapor (mm) 

RH mean relative humidity (%) 

RHmax Maximum relative humidity (%) 

S monthly average daily bright sunshine duration (h) 


So monthly average maximum possible daily sunshine 
duration (h) 

Suh sunshine duration that take into account the natural 
horizon of the site (h) 

ST mean soil temperature (°C) 

T mean air temperatures (°C) 

Tk mean air temperature (K) 

T wax mean maximum daily temperatures (°C) 

Tin mean minimum daily temperatures (°C) 

AT Temperature difference (°C) 

ATm monthly mean AT (*C) 

ta mean dew point temperature (°C) 

V mean water vapor pressure (hpa) 

w atmospheric precipitable water vapor per unit 
volume of air (cm) 

Z altitude of site (km) 

es(T) saturation vapor pressure at the air temperature 
T (kPa) 


Greek letters 


ò solar declination (°) 

0 solar zenith angle (°) 

p latitude of the location (°) 

L latitude of the location (rad) 
A longitude of the location (°) 


Ws sunset hour angle (°) 

hy the noon solar altitude angle of the sun (°) 

T atmospheric transmittance 

Tt, max maximum (cloud-free) daily total transmittance at a 
location 

Tf,max proportion Of Ttmax Observed on a given day (cloud 
correction) 


of average daily global radiation to the corresponding value on a 

completely clear day (H/H.) at a given location and the ratio of 

average daily sunshine duration to the maximum possible sun- 

shine duration. 

H S 

m=’ +b (=) (1) 
A basic difficulty with Eq. (1) lies in the definition of the term 

Hc. Prescott [4] and the others have modified the method to base 

it on extraterrestrial radiation on a horizontal surface rather than 

on clear day radiation and therefore proposed the following 

relation: 

H S 

Ho = a+b (=) (2) 
The values of the monthly average daily extraterrestrial 

irradiation (Ho) can be calculated from the following equation [5]: 


Wh 24 360n 
Ho (45) = pu [140.033 cos( 365 ) 


[cosp COSÓ sinc)», + €, sinp sinà| (3) 


T 
180 

The solar declination (6) and the mean sunrise hour angle (c) 
can be calculated by following equations [5]: 


" . [360 
6 — 23.45 sin [ses (284+ »| (4) 


ws = cos! (- tanq tanó) (5) 


For a given month, the maximum possible sunshine duration 
(monthly average day length, Sy) can be computed by using the 
following equation [5]: 


So = 15 Ds (6) 

Although the Angstróm-Prescott [4] equation can be improved to 
produce more accurate results, it is used as such for many applica- 
tions. Some of the regression models based on the Angstróm-Prescott 
model which proposed in literature are given as follows: 


- Jain model for Italy [6] 


H 
Ho 


— 0.1774 o.eo2() (7) 
So 


- El-Metwally model for Egypt [7] 


H S 
HU 0228-0527 (2) (8) 


- Bakirci model for Turkey [8] 
H S 
m 0.2786+ 0.41 60(,) (9) 


- Alsaad model for Amman, Jordan [9] 


H S 
Hs = 0.1744 0.615(=-) (10) 
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Jain and Jain model for Zambian [10] 


H 


Bia =0.240+0.513( $) (11) 
Ho 


So 


Katiyar et al. model for India [1] 


H S 
Ho 0.2281 4- 0.5093 (=) (12) 


Lewis model for state of Tennessee, U.S.A. [11] 


A — 0.144 057(=) (13) 

Tiris et al. model for Turkey [12] 

H 0484 0.62(=) (14) 
0 So 

Almorox and Hontoria model for Spain [13] 

A — 0.21704- 0.5453 (=) (15) 


Raja and Twidell model for Pakistan [14,15] 


H 
Ho 


= 0.3354 0367(,) (16) 
So 


Li et al. model for Tibet, China [16] 


H S 
Ho = 0.2223+0.6529 (s) (17) 


Said et al. model for Tripoli, Libya [17] 


H 
Ho 


— 0.2154 0527(=) (18) 
So 


Ulgen and Ozbalta model for Izmir-Bornova, Turkey [18] 


H S 
Ho -02424.«05014(2) (19) 


El-Sebaii et al. model For Jeddah, Saudi Arabia [19] 


H S 
Hy =7281+3.78 (=) (20) 


Rensheng et al. model for 86 stations in China [20] 


H S 
Ho =0.176+0.563 (=) (21) 
Ahmad and Ulfat model for Karachi, Pakistan [21] 


H 
Ho 


= 0.324+0.405 ( =) (22) 
So 


Jin et al. model for 69 stations in China [22] 


H S 
Hi 0.1332+0.6471 (5) (23) 


Akpabio and Etuk model for Onne region (within the rainforest 
climatic zone of southern Nigeria) [23] 


H S 
HU 0.23 +0.38 (=) (24) 


- Ulgen and Hepbasli model for Ankara, Istanbul and Izmir in 


Turkey [24] 


H S 
Ho = 02671404754 (=) (25) 


Aras et al. model for twelve provinces in the Central Anatolia 
Region of Turkey [25] 

H 
Ho 


= 0.3078+0.4166 (s) (26) 
0 

Togrul and Togrul model for Ankara, Antalya, Izmir, Yenihisar 

(Aydin), Yumurtalik (Adana) and Elazig in Turkey [26] 


H 
Ho 


= 0.3184 0.449( $) (27) 
So 


Kholagi et al. model for three different stations in Yemen [27] 


H S 
Ho = 019140571 (=) (28a) 
H S 
Ho = 0297+0.432 (=) (28b) 
H S 
Ho = 0262-0454 (=) (280) 


Katiyar et al. model for Jodhpur, Calcutta, Bombay, Pune in 
India, respectively [1] 


H s 
p; = 02276-05105 (=) (29a) 
H ..02623403952( 2 (29b) 
Ho s 
H s 
p; = 02229405123 (=) (290) 
H ..022864.0.5309 (s) (29d) 
Ho So 


- Jain model for (Salisbury, Bulawayo and Macerata, Italy), 


respectively [28] 


H S 

Hg — 0.3134 o474($) (30a) 

H S 

B = 0.307 -- 0.488 (=) (30b) 

H = 0.3094 0.599(=-) (30c) 
0 So 


Veeran and Kumar model for two tropical locations Madras 
and Kodaikanal in India, respectively [29] 


H 
0 


H S 
Ho =0.27+0.65 (=) 


= 0.344 0.32(=) (31a) 
So 


(31b) 


Chegaar and Chibani model for Algiers and Oran, Beni Abbas 
and Tamanrasset, Algeria [30] 


H 030910368 (=) 
Ho 


5 (32a) 
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H S 
p; = 0367 +0.367 (=) (32b) 
H S 

Ho 7023340591 (=) (32c) 


- Ampratwum and Dorvlo model for Seeb and Sallalah weather 
stations in Oman, respectively [31] 


H = 0.33264 03110(=) (33a) 
0 So 

H S 

Hs = 0.2418 4- 0.3555 (=) (33b) 


2.1.2. Model 2: Glower and McCulloch model 

Glower and McCulloch proposed the following equation, 
which takes into account the effect of latitude of the site, @, as 
an additional input and is valid for ( < 607)[32]: 


H S 
Ho =acosp+b (=) (34a) 
H S 
Ho = 0.29 cosp+0.52 (s) (34b) 
Some other corresponding modified relations are: 
- Ulgen and Hepbasli model for Izmir, Turkey [33] 
H S 
— = 0.3092 coso 4- 0.4931 | — (35) 
Ho So 
- Raja and Twidell model for Pakistan [14,15] 
H = 0.388 cosp+ 0.367 2 (36a) 
Ho So 
H _ 0.388 cosp + 0.407 (=) (36b) 
Ho So 


2.1.3. Model 3: Page model 

Page has provided the coefficients of the modified Angstróm- 
type model [4], which is claimed to be applicable anywhere in the 
world [34]: 
H 
Ho 


=0.234 o4s(2 ) (37) 
So 


2.1.4, Model 4: Rietveld model 

Rietveld by using measured data collected from 42 stations 
located in different countries, has proposed a unified correlation 
to compute the horizontal global solar radiation. Rietveld's model, 
which is claimed to be applicable anywhere in the world, is given 
in following equation [35]: 
H 
Ho 


—0.184 oe2(2) 38) 
So 


Rietveld also examined several published values of a and b 
coefficients of the Angstróm-Prescott model [4] and noted that 
constants a and b are related linearly to the appropriate mean 
value of S/Sp as follows [35]: 


a=0.10+0.24 (=) (39a) 
So 
S 

b=0.38+0.08 (=) (39b) 
0 


2.1.5. Model 5: Dogniaux and Lemoine model 

Dogniaux and Lemoine proposed the following correlation, 
where the regression coefficients of the Angstróm-Prescott model 
[4] seem to be as a function of the latitude of the site [36]: 


a = 0.370224 (0.00506 (=) -0.00313 p+0.32029 (=) (40) 
0 


0 0 
They also in same year obtained the specific monthly correla- 
tions which are listed in Eq (41a)-(411) [36]: 


January B = (—0.00301 + 0.34507) + (0.00495 + 0.34572) (=) 
0 0 


(41a) 


February 
0 0 


(41b) 


= = (-0.00255p + 0.33459) + (0.004579 + 0.35533) (= 


March + = (—-0.00303— +0.36690) + (0.00466¢ + 0.36377) (5) 
0 


(41c) 

~ H S 
April — = (-0.00334p+0.38557) + (0.004569 + 0.35802) (=) 
0 
(41d) 


May -=(-0.00245p+0.35057) + (0.00485 + 0.33550) (=| 


June = = (—0.00327 + 0.39890) + (0.00578 +.0.27292) 8 
0 


(411) 


July 7 = (-0.00369¢ +0.41234) + (0.00568 +0.27004) (=) 
0 0 


(41g) 


August a = (—0.00269¢ + 0.36243) + (0.00412 +0.33162) (=) 
0 0 


(41h) 


September a = (—0.00338¢ + 0.39467) + (0.00564 + 0.27125) (=) 
0 o 


(41i) 


October = = (—0.00317 4-0.36213) + (0.00504 4-0.31790) (s) 


0 0 


(41j) 


November E = (—0.00350¢ + 0.36680) + (0.005239 + 0.31467) (=) 


0 0 


(41k) 


December ED = (—0.00350¢ + 0.36262) + (0.005590 + 0.30675) (= 
0 o 


(411) 


2.1.6. Model 6: Kilic and Ozturk model 
Kilic and Ozturk calculated the a and b regression coefficients 
of Angstróm model [4] for Turkey [37]: 


a = 0.103 + 0.000017 Z - 0.198cos(q—ó) (42a) 


b = 0.533—0.165 cos(q—ó) (42b) 
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2.1.7, Model 7: Benson et al. model 

Benson et al. proposed two different formulations for two 
intervals of a year depending on the climatic parameters to 
estimate the global solar radiation [38]: 

For January-March and October-December 


H S 

Ho =0.18+06(=) (43a) 
For April-September 

H S 

m 0.24 4- 0.53 (=) (43b) 


2.1.8. Model 8: Ogelman et al. model 

Ogelman et al. expressed the ratio of global to extraterrestrial 
radiation by a second order polynomial function of the ratio of 
sunshine duration [39]: 


H S sy? 
Ho (=) l «(s) P 
H 049540676( 2-) -o142( 2 Í (44b) 


Many authors in all over the world applied this model and 
determined the regression coefficient of this model for particular 
location of interest as follows: 


- Akinoglu and Ecevit model for Turkey [40] 


2 
H — 0.1454 0.845(=) 0.280( $) (45) 
0 So So 
- Almorox and Hontoria model for Spain [13] 
H S sy? 
Hm 0.1840 -- 0.6792 (&) —0.1228 (=) (46) 
- Bakirci model for Turkey [8] 
H S SA 
m 0.2545 4- 0.5121 (=) —0.0864 (=) (47) 


- Tasdemiroglu and Sever model for Ankara, Antalya, Diyarbakir, 
Gebze, Izmir and Samsun in Turkey [41] 


H 
Ho 


S sy? 
=0.225+0.014 +0.001 (2 (48) 
So So 


- Yildiz and Oz model for Turkey [42] 
H 


2 
H7 0.2038 + 0.9236 (s) —0.3911 (2) (49) 
0 


So So 
- Aksoy model for Ankara, Antalya, Samsun, Konya, Urfa and 
Izmir, Turkey [43] 


H 
0 


S sy? 
= 0.148 +0.668 0.079 (50) 
So So 


- Said et al. model for Tripoli, Libya [17] 


H S sy? 
H; = 0.1 + 0.874 (=) —0.255 (s) (51) 


- Ulgen and Ozbalta model for Izmir-Bornova, Turkey [18] 


H S sy? 
qm 0.0959+0.9958 (=) —0.3922 (=-) (52) 


- Togrul and Togrul model for Ankara, Antalya, Izmir, Yenihisar 
(Aydin), Yumurtalik (Adana) and Elazig in Turkey [26] 


H S sy? 
Hy =0.154141.1714( =) -0.705 (=) (53) 


- Ahmad and Ulfat model for Karachi, Pakistan [21] 


H S sy? 
HU 0348032027) -- 0.070 (=) (54) 


- Tahran and Sari model for Central Black Sea Region of Turkey [44] 


H S sy? 
HU 0.1874 4- 0.8592 (=) —0.4764 (s) (55) 


- Jin et al. model for 69 stations in China [22] 


H 


2 
i, = 0.1404+ 0.6126 (=) +0.0351 (=) (56) 
0 


So So 


- Aras et al. model for twelve provinces in the Central Anatolia 
Region of Turkey [25] 


H 
Ho 


2 
= 0.3398 + 0.2868 (=) +0.1187 (=) (57) 
So So 


- Ampratwum and Dorvlo model for Seeb and Sallalah weather 
stations in Oman, respectively [31] 


H S sy? 
Ho = 0.9428— 1.2027 (=) H 0.9336(%-) (58a) 
H S sy? 
m7 0.1971 +0.6297 (=) —0.2637 (=) (58b) 
2.1.9. Model 9: Bahel et al. model 
Bahel et al. reported the following relationship [45]: 
H S 
He =0.175+0.552 (=) (59) 


2.1.10. Model 10: Zabara model 

Zabara correlated monthly a and b values of the Angstróm- 
Prescott model [4] with monthly relative sunshine duration (S/So) 
as a third order function and expressed the a and b coefficients as 
[46]: 


S sy? sy? 
a=0.395—1.247( —)+2.680( — | —1.674( — (60a) 
So So So 

S sy? sy? 
b=0.395+1.384 3.249 22.055 (60b) 
So So So 


2.1.11. Model 11: Bahel model 

Bahel developed a worldwide correlation based on bright 
sunshine hours and global radiation data of 48 stations around 
the world, with varied meteorological conditions and a wide 
distribution of geographic locations [47]. 


H S sy? sy? 
Hp =016+087(3) oe1(2) -034( =) (61) 


2.1.12. Model 12: Gopinathan model 
Gopinathan suggested a and b regression coefficients of 
Angstróm-Prescott model [4] as a function of elevation (Z) and 
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sunshine ratio (S/So) for estimation of the global solar radiation [48]: 


a=0.265+0.07Z— 0.135 (s) (62a) 
0 
b=0.401—0.108Z+0.325 (=) (62b) 
0 


Gopinathan also reported the following correlations [49]: 


A = |—0.309+ 0.539 cosq—0.0693Z + 0.290 Ed 
Ho So 


+ |1.527-1.027 coso 4- 0.0926Z o35e($)] (=) (63) 
0 0 


2.1.13. Model 13: Newland model 
A linear-logarithmic model, has been used by Newland to 
obtain the best correlation between H/Hoand S/Sy [50]: 


H S S 

m =e+o(3 -) + clog ( $ ) (64a) 

H S S 

7, - 034040 (5) +0.17108 (=) (64b) 
Some applications of this model are: 

- Bakirci model for Turkey [8] 
H = 0.3925 - 0.2877 2 +0.1527 log 2 (65) 
Ho So So 


- Ampratwum and Dorvlo model for Seeb and Sallalah weather 
stations in Oman, respectively [31] 


H S S 

m^ —1.1931-- 1.8641 (=)- —1.2544 log (=) (66a) 
Aa = 0.4337 + 0.1430 : 4- 0.0861 log > (66b) 
Ho So So 


2.1.14. Model 14: Monthly specific Rietveld model (Soler model) 
Soler applied Rietveld’s model to 100 European stations and 
gave the following specific monthly correlations [51]: 


January a = 0.184 0.56(,) (67a) 
February in = 0.20-- 0.60 (=) (67b) 
0 
March d = 0.22 +0.58 s (67c) 
Ho So 
April zi = 0.20-- 0.62 ( =) (67d) 
0 So 
H S 
May Ho = 0.244 0.52 (=) (67e) 
H S 
June m 0.244- 0.53 (=) (67f) 
H S 
July m = 0.23 4- 0.53 (s) (678) 
August H = 0.22+ 0.55 s (67h) 
Ho So 


September H = 0.204- 0.59 S (671) 
Ho So 

October H — 0.194-0.60 3 (67j) 

Ho So 

November H = 0.174 0.66(=) (67k) 
0 So 
H S 

December — = 0.18 +0.65 (671) 
Ho So 


Also, they reported the regression coefficients of a and b as 
follows [51]: 


a=0.179+-0.099(=-) (68a) 
0 


2 
b=0.1640+0.1786(5,)- 1 ooss (3) (68b) 


0 So 


2.1.15. Model 15: Luhanga and Andringa model 
Luhanga and Andringa derived their own model as follow [52]: 


H 
0 


= 0.241 + 0.488 ( S ) (69) 
So 


2.1.16. Model 16: Louche et al. model 
Louche et al. presented the model below to predict global solar 
radiation [53]: 


H 


ji, =0.206+0.546( =) (70) 
Ho 


So 
Furthermore, the Angstróm-Prescott model [4] has been 
modified through the use of the ratio of (S/S,.)instead of (S/So) 
by Louche et al. model which is presented as follows [53]: 


H S 

= a+b(3.) 71a 
Ho Shh ( ) 
{1 _ 0-8706 0.0003 (71b) 
Sin So 


2.1.17. Model 17: Samuel model 
Samuel has correlated H/Howith S/Sg in the form of a third 
order polynomial equation [54]: 


H S sy? sy? 
Hy (5) i (S) +a(3,) 


H S sy? sy? 
Hp =-0:14+252(5,) 3.71(5,) -224($) 


Some investigator used this model to estimate global solar 
radiation in different locations all over the world. Some of these 
models are: 


(72a) 


(72b) 


- Ertekin and yaldiz model for Antalya, Turkey [55] 
H S sy? sy? 
H^ —2.42754- 11.946 (=)- —16.745 (=) +7. sz (s) 
(73) 
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- Almorox and Hontoria model for Spain [13] 


H S sy? sy? 
Hp =0230+03809 (7) +0.4694( $) 0.3657( =) (74) 


- Lewis model for state of Tennessee, U.S.A. [11] 


H S ae gx 

Hi 0.81—3.34 (=) +7.38 (=) —4.51 (=) (75) 
- Ulgen and Hepbasli model for Izmir, Turkey [33] 

H S sy? S3 

Hg = 0.2408 + 0.3625 (s) - 0.4597 (s) 0.3708 (=) (76) 


- Togrul and Togrul model for Ankara, Antalya, Izmir, Yenihisar 
(Aydin), Yumurtalik (Adana) and Elazig in Turkey [26] 


2 3 
H = 0.17964 09813.) 0.2058 (7) 0.2657(.) (77) 


Ho So So So 


- Ulgen and Hepbasli model for Ankara, Istanbul and Izmir in 
Turkey [24] 


2 3 
ii = 0.2854-- 0.2591 (s) + 0.6171 (s) 0.4834(5,) (78) 


Ho So So So 


- Tahran and Sari model for Central Black Sea Region of 


Turkey [44] 
H S sy? 5\3 
Ho = 0.15204 11334(2.) 1.1126(=) | 0.4516( =) (79) 


- Jin et al. model for 69 stations in China [22] 


H s sy? S43 
Ho =0.1275+0.7251 (=) 0.2299( =) | 0.1837(7,) (80) 


- Aras et al. model for twelve provinces in the Central 
Anatolia Region of Turkey [25] 


2 3 
H _0.4832-06161(2)+1.8932(2) -1.0975( 
Ho So So 


- Rensheng et al. model for 86 stations in China [20] 


H S sy? E 
H =0.150+1.145 (=) EX ) -0.963(=) (82) 


0 0 So 


- Bakirci Model for Erzurum, Turkey [56] 


H S sy? S 
Hi — 0.6307 0.7251( =) +1.2089( $) 0.4633 ( =) (83) 


2.1.18. Model 18: Gopinathan and Soler model 
Gopinathan and Soler suggested following linear equations for 
locations with latitudes between 60°N and 70°N [57]: 


H S 
p; = 01538-07874 (=) (84a) 
H S 
Hy -01961 +07212( Š) (84b) 


2.1.19. Model 19: Raja model 
Based on Bennett's formula [58], Raja proposed the following 
insolation-sunshine relation [59]: 


H Z 
Ho (0.368-0.125 1-z)] 


+ [0567-0018 (1-7) —0.211 coso 2 (85a) 


Za So4 


where Za = 8000 m and So4 is the 4° corrected day length used to 
compensate the finite threshold of the Champbell-Stokes sunshine 
recorder and is calculated as follows: 


2 1 [(sin4^—singsinó) 
=j | cosócosq | 


Soa (85b) 


2.1.20. Model 20: Coppolino model 
Coppolino developed a power function and incorporated a 
trigonometric term to estimate global solar radiation [60]: 


(a) Power model 
b 
ño e (s) as 


- Ampratwum and Dorvlo model for Seeb and Sallalah 
weather stations in Oman, respectively [31] 


0.4253 
5 = g- 04470 (=) (87a) 
0.3588 
ie — g 05561 (=) (87b) 
(b) Power-trigonometric model 
b 
M =e! (=) (sinh;) (88) 


The noon solar altitude angle of the Sun (h,) is the comple- 
ment of zenith angle and can be calculated from [60]: 


- Coppolino employed this model to estimate global solar 


radiation in 34 Italian stations. The regression coefficients 
of this model are presented as follows [60] 


H S 945 ; 0.05 S 
Hs — 0.67 (=) (sinh) 0.15 < (=) < 0.90 (90) 


- Ampratwum and Dorvlo model for Seeb and Sallalah 
weather stations in Oman, respectively [31] 


0.5239 
i =e 04135 (=) (sinh, 29 (91a) 
0.3571 
M — g-05598 (=) (sinh) 00296 (91b) 
2.1.21. Model 21: Tiris et al. model 
Tiris et al. reported their own correlations as follow [61]: 
& — 0.22624 0.418 (<) (92) 


2.1.22. Model 22: Ampratwum and Dorvlo model 
Ampratwum and Dorvlo have suggested a logarithmic type 
model as [31]: 


H S 
H; =a+blog(3-) (93) 
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Some other developed logarithmic correlations are: 


- Almorox and Hontoria model for Spain [13] 


A 0.6902 + 0.6142 log n (94) 
Ho So 


- Ampratwum and Dorvlo model for Seeb and Sallalah weather 
stations in Oman, respectively [31] 


H 0.6376+0.2490 log( > (95a) 
Ho So 
H _0.5612+0.1412 log( > (95b) 
Ho So 

- Bakirci model for Turkey [8] 
2T = 0.6446 + 0.4842 log 2 (96) 
Ho So 


- Togrul and Togrul model for Ankara, Antalya, Izmir, Yenihisar 
(Aydin), Yumurtalik (Adana) and Elazig in Turkey [26] 


H 0.698+0.2022 log(> (97) 
Ho So 


2.1.23. Model 23: Klabzuba et al. model 
Klabzuba et al. developed the following relationship between 
the measured H and the relative sunshine duration [62]: 


S (S 2 
H=a4 b(=) «cla (2)]e e) (98) 


- Klabzuba et al. model for Hradec Králové solar observatory, 
Czech Republic [62] 


H=7.19+0.258 (=) 9.28 x 10-9 [8 4229 |m-1747* (99) 
0 0 


2.1.24. Model 24: Togrul et al. model 

Togrul et al. applied the various regression analyses to inves- 
tigate the relation between monthly mean S/So ratio and con- 
stants a and b in Angstróm model [4] as follow [63]: 


E 5 
a- is): b=k+1(5) (100a) 
a=itn(5,) Hj; b- kin (s )+m (100b) 
(Ss EXT S sy? 
a=i+j(5,) Ze : b=t+m(,) (s) (100c) 


eco) als) (S) roms) (o) (5) 


(100d) 
a=i+j $ Ek sy l sy m si 
=i+j z) * s + s + 5 
S sy? sy? sy? 
ld (5) '(z) ' (5) | o(5) 


- Togrul et al. obtained the empirical coefficients of their 
models for Elazig, Turkey as follows [63] 


(100e) 


a= 0.2816 0.1049( E 
So 


|: b = 0.4762 0.0174( =.) 
So 


(101a) 


a=-—0.0344In (=) +0.1982; b= —0.0201In B +0.4562 
0 0 


(101b) 


b — 0.6248—0.7033 


a = 0.1950-- 0.2943 (s) —0.3706 
0 


2 
5) (1010) 


So 


0.9913 


a=0.2646-0.2958 ( =) -0.9671 (=) 08749 B? 
( 


b= 0.5403 +-0.0149( 3 
So 


(101d) 


S S^ 
a — 0.22264-0.2841 ( — | -1.3653( — 
So So 


SA? sy? 
+2.6661(¿,) -1.7997(,) 


S sy? 
b=0.561-0.2715 -0.1606[ > 
So So 


sy? sy? 
-0,6839( =) +0.8888(5-) 


(101e) 


2.1.25. Model 25: Elagib and Mansell model 

Elagib and Mansell developed new techniques for predicting 
solar radiation based on sunshine hours and geographical para- 
meters [64]: 


Lb) v. 
j, =a+bl+cz+d( $) (102b) 
i =a+bl+e(=) (102c) 
E =a+b2+c(=) (102d) 


Some researchers calibrated the Elagib and Mansell models for 
different locations. Some of these models are as follows: 


- Togrul and Togrul model for Ankara, Antalya, Izmir, Yenihisar 
(Aydin), Yumurtalik (Adana) and Elazig in Turkey [26] 


H _ 9.3396 exp|0.8985 (> (103) 
Ho So 

- Jin et al. model for 69 stations in China [22] 
H = 0.0855 4- 0.0021 4-0.03Z + 0.5654 ES (104) 
Ho So 

- Rensheng et al. model for 86 stations in China [20] 
a = 0.1224 0.001L 42.57 x 10°7Z+0.543 2 (105) 
Ho So 


2.1.26. Model 26: Monthly specific Elagib and Mansell model 
Elagib and Mansell have investigated the possibility of estab- 

lishing monthly-specific equations for estimating global solar 

radiation across Sudan. The best performing equations for each 
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month are given as following [64]: 


January ue = 0.1357 +0.3204L + 0.0422Z 4- 0.4947 (=) (106a) 
0 


0 


February = = 0.1563 +0.3166L+0.1006Z + 0.4593 (=) 
0 


So 
(106b) 
H g \ 07263 
March He ome) (106c) 
4, H S 
April = 0.1640 + 0.0397Z + 0.5773 (106d) 
Ho So 
H S sy? 
May — =0.0709+0.8967| —]-0.2258| — (106e) 
Ho So So 
2 
June H _ 00348415078 2 —0.8246 2 (106f) 
Ho So So 
H S 
July H = 0.3205 + 0.1444L + 0.0782Z + 0.2916 [as (106g) 
0 0 
H S 
August H, =0-2720+0.0369L4-0.101724+0.3888 $ (106h) 
0 0 
H S sy? : 
September =~ — —0.3710-4- 2.57831 — | -1.6788| = (106i) 
Ho So So 
H S : 
October H = 0.1593-—0.1043L + 0.0609Z + 0.5916 5 (106j) 
0 0 
H S 
November = 0.1786 +0.0199Z + 0.5441 | — (106k) 
Ho So 
H S 
December H. = 0.1714 -- 0.1329L + 0.0482Z -- 0.5015 s 
0 0 
(1061) 


2.1.27. Model 27: Togrul et al. model 

Togrul et al. developed some statistical relations to estimate 
monthly mean daily global solar radiation by using clear sky 
radiation. The correlations have been developed by employing 
both ratios of (S/So) and (S/Syn)[65]. 

where Spn is the monthly mean sunshine duration that take 
into account the natural horizon of the site and can be calculated 
by following equation: 


" = Bd 1.0.0003 (107) 
ni 
For summer 

2 
H 05771406843 (s) —0.3544 (=) (108a) 
C 0 0 
H  043564.13939( 3) 14988 (3 +05054(2) cosb 
MM C RD RM A 
= — 0.7295exp [02633 (&)] (1080) 
Cc 

0.1688 

= E 09337 (5) (108d) 
£ 
H 0577407825( 2-)-o46x1( 5) (108€) 
He i = Snh : Sun 


2 3 
a = 0.4356 + 1.5935 (=) —1.9584 (=) +0.8893 (=) 


He Snn Sun Suh 

(108f) 
H 5 \ 0.1689 
— —0.9552| — 108 
Hc (=) ( 8) 

For winter 

H S sy? 
m^ 0.2384 4- 2.1498 (=) —1.8664 (=) (109a) 


H S sy? sy? 
p, =03081+1.4933( Š) 0.0839(3-) 1.4513 (>) (109b) 


H 05451 exp (0.7981 8] a099 
c So 

H S 0.3204 

H; = 10302 (=) (109d) 

H 0.2385+2.4604 (=) -2.445 (2) (109e) 

Hc EL f Snh f Sin 


H S sy? sy? 
— =0.3077+1.7137 (=) —0.1251 (5) —2.1612 (=) 
Hc Sah Sun Sun 
(109f) 


0.3202 
H 10755 (=) (109g) 


2.1.28. Model 28: Almorox and Hontoria model 
Almorox and Hontoria proposed the following exponential corre- 
lation to estimate global solar radiation from sunshine hours [13]: 


H S 
g t5 exp( =) (110) 
- Almorox and Hontoria model for Spain [13] 
H _0.0271+0.3096 exp( > (111) 
Ho So 
- Bakirci model for Turkey [8] 
H _ 9.0963+0.2337 exp( 3. (112) 
Ho So 


2.1.29. Model 29: Jin et al. model 

Jin et al. proposed the following models by using solar 
radiation data and some geographical parameters like latitude 
and altitude [22]: 


H 


S 
Ho — a4- bcosq 4- cZ4 a(=) (113a) 
H S 
ciis =(a+bL+c2)+d+eL+12)(5,) (113b) 
Ho So 
H S 
ID (a-- bcosq +cZ) + (d-- ecoso -- fZ) (s) (1130) 

2 

H — (a -- bL-- cZ) - (d 4- eL a(S) H(g+hL4 D$) (113d) 
Ho So So 


H S 
m (a+ bcosq +cZ) + (d-- ecosq -- fZ) (s) 
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2 
+ (8+hcosp+1iZ) (=) (113e) 


Some implications of these models are: 


- Jin et al. model for 69 stations in China [22] 


H 


n 2.1186—2.0014 cosp 4- 0.0304Z + 0.5622 (=) (114a) 
0 


0 


H 
Ho 


= (0.1094+0.0014L+ 0.02122) 


+ (0.5176 -+0.0012L + 0.01502) (=) (114b) 
0 


H _ (1.8790-1.7516 cos +0.0205Z) 


+ (1.0819—0.5409 cos + 0.01692) (=) (114c) 


0 


H = (0.0218 +0.0033L+ 0.04432) 
0 


+ (0.9979—0.0092L 0.08522)( =) 
0 


2 
+(-0.5579+0.0120p +0.1005Z) (=) (114d) 
0 


—— = (4.2510-4.1878 cosp+0. 7Z 
H 4.2510—4.1878 0.043 
0 
S 


+ (710.5774 - 11.4512 cosp—0.0832Z) (5) (114e) 
0 


- Rensheng et al. model for 86 stations in China [20] 


H 0280-0141 cosq 4-2.60 x 107Z--0.542 (=) (115a) 
Ho So 
H _ (0.:329-0.196 2.2 x 10°7Z 
Ho = ( ) —0. cosp+2.2 x ) 
$ (0.457 --0.097 cosp +6.72 x 10?z) (=) (115b) 
0 


2.1.30. Model 30: El-Metwally model 
El-Metwally proposed a non-linear correlation between clear 
index (H/Ho) and relative sunshine (S/S) as follows [7]: 


H (^ ¿(us/5) 


Ho = a (116) 
- El-Metwally model for Egypt [7] 
= — 0.713(/5/5)) (117) 


0 


2.1.31. Model 31: Almorox et al. model 

Almorox et al. reported the monthly-specific equations for 
estimating global solar radiation from sunshine hours for Toledo, 
Spain as given by following correlations [66]: 


S 


January Es = 0.285 + 0.444 (5) (118a) 
Ho So 
H S 
February — = 0.272 + 0.465 | — (118b) 
Ho So 
H S 
March | —— — 0.291 4-0.491| = (118c) 
Ho So 


April = = 0.266+0.495 ( 3 ) (118d) 
Ho So 
H S 
May 1 =0.286+0.475 (=) (118e) 
Ho So 
June X = 0.311-- 0.439 (=) (118f) 
Ho So 
H S 
July — = 0.329 +0.406 (s) (118g) 
Ho So 
H S 
August — 0.3134 0.410( ) (118h) 
Ho So 
H S ; 
September = 0.271 -- 0.479 (118i) 
Ho So 
H S ] 
October — —0.259-- 0.465 | — (118j) 
Ho So 
H S 
November —=0.279+0.431| — (118k) 
Ho So 
H S 
December —=0.282+0.428| = (1181) 
Ho So 


2.1.32. Model 32: Rensheng et al. model 

Rensheng et al. have suggested new equations, based on the 
Angstróm model [4] to estimate global radiation. These equations 
are as follows [20]: 


H S sy? sy? 
Hs — (a4- bL-- cZ)4 a(s,) (s) f(s) (119a) 
H S sy? sy? 
H^ (a+ bcosq4 2) +d(,) +e(,) f(s) (119b) 
H i i i S 
He = (a+ bcosp+cZ) + (d+ecosp+fZ) s 
sy? sy? 
+ (8+hcosq +iZ) (=) + (¡+kcosp+IZ) e (1190) 
0 0 
2E =(a+bL+cz+di)+e(3,) (119d) 
Ho So 
d ec b Z4-dÀ 3 119 
pg; "ut COSQ -- CZ - dA) +e 5 (119e) 
H S 
=(a+bL+cZ+d2)+(e+fL+8Z+ n»( ) (119f) 
Ho So 
H S 
Hy = (a+ bcosq -- cZ 4- d4) + (e--fcosq - gZ - hA) (=) (119g) 


H _ 2 A SN 
He = (a+bL+cZ+di+ed )+#( =) (s) ' n(5,) (119h) 


H " S sy? sy? 
HU (a+beosp-+cZ+di+ed )+* (5) +8(5,) +h(3,) 
(119i) 


—= (a+bcosp+cZ+d1+e2?) 


+ (f +8cosp+hZ+11+ i?) (=) 
So 


810 


2 
+ (k +Icosp+mZ+n2+ 07?) (=) 


3 
E (p+ qcosQ--TZ -sÀ4- ti”) (=) 
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So 


So (119) 


- Rensheng et al. employed their models to obtain following 


correlations using data from 1994 to 1998 at 86 stations in 
China [20]. 
H (0.109 +0.0011+2.41 x 1072) +1 o29( 5. 
Ho MX” mE B i So 
sy? sy? 
1.216 +0.787 (120a) 
So So 
Hz (0.234-0.112 cosp+2.43 x 1077z) 
Ho 
S sy? sy 
+1.026 (=) 1.209 (=) + 0.782 (=) (120b) 
LL (0.336-0.233 cosp+2.64 x 107Z) 
Ho 
+ (—0.670+2.140 cosp—0.1Z) (=) 
0 
S 2 
+ (2.744—5 cosq -- 0.3Z) (=) 
0 
S 3 
+ (-1.638+3.042 cosp—0.2Z) (5) (1200) 
0 
LL (0.117-0.0011+2.59 x 1022-411 x 10-57) 
Ho 
+0.543 (s) (120d) 
So 
f (0275-0141 cosp+2.63 x 1077Z 
0 
= S 
+4.27 x 10 i) 4-0.542 ( > (120e) 
So 
- = (0.094 +0.002L+ 2.27 x 107?Z+0.00014) 
0 
+ (0.586—0.0008L.- 5.36 x 10 *Z—0.00022) (=) 
0 
(120f) 
Lm (0.313-0.195 cosp+2.28 x 10-7Z+0.00012) 
0 
+ (0.476+0.097 cosp +5.69 x 107°Z—0.00024) (=) 
0 
(120g) 
ui - (0.370 +0.0007L+2.44 x 10?Z—0.0054 4-2.24 x 10577) 
0 
-1026( 5) 1208( S 40783 (5 (120h) 
E (2) i (s) E (s) 
E - (0.426-0.087 cosp+2.44 x 107?Z—0.0044+ 1.86 x 10757) 
0 
¿1024 (12045 4077905 (120i) 
E So : So) * So 
ue = (1.012-0.141 cosp+2.84 x 10-?2—0.0144 
0 


462 x 10772) + (-4.061 +1.740 cos 


—0.1Z+ 0.0694—0.0034?) (s) + (12.402—3.867 coso 
0 


2 
+0.32-0.1997.+-0.00097?) (=) + (-9.442 
0 


3 
+2.115 cosp—0.2Z+ 0.164/—0.00087?) (=) 
0 


(120j) 


2.1.33. Model 33: Sen model 

Sen proposed a nonlinear model for the estimation of global 
solar radiation from available sunshine duration data. This model 
is an Angstróm type model with a third parameter appears as the 
power of the sunshine duration ratio [67]: 


(121) 


— El-Sebaii et al. calibrated this model For Jeddah, Saudi Arabia [19] 


2.344 
H __0.864+1.862 (=) 


Es = (122) 


2.1.34. Model 34: Bakirci model 

Bakirci reported the original Angstróm-type equations including 
the linear, second-order and fifth-order polynomial relationships 
between the monthly average values of (H/H.) and (S/So) as 
follows [68]: 


H S 
H, = 07836 0.0460( 7.) (123a) 
H _4.9192-1.0547(*) +0.9661 (2 : 123b) 
He So)" ` So ( 
H __11.2254+128.010(2)-516.900(2 * +994.730 sy 
He 07 So So f So 
sy? sy? 
-920.350(2) +329.93( > (1230) 
So So 


2.1.35. Model 35: Bakirci model 
Bakirci developed following models which is a derivation of 
the modified Angstróm-type equation [8]: 


a) Linear exponentia 

(a) Li p ial 
g a+b Na cex 3 
Ho |  VMSo ` p So 


Some other relations are: 


124) 


- Muzathik et al. model for state of Terengganu, Malaysia [69] 


H = 0.19490 + 0.4771 = +0.02994exp = (125) 
Ho So So 

- Bakirci model for Turkey [8] 
2t = 0.3448 + 0.5636 2 —0.0838exp = (126) 
Ho So So 
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(b) Exponent 
H sy? 
H^ a (s) (127) 
- Bakirci model for Turkey [8] 
H 5 04333 
H^ 0.6660 (s) (128) 


- Togrul and Togrul model for Ankara, Antalya, Izmir, Yenihisar 
(Aydin), Yumurtalik (Adana) and Elazig in Turkey [26] 


0.4146 
He =97316(g,) 


HU a (129) 


2.2. Cloud-based models 


Clouds and their accompanying weather patterns are among the 
most important atmospheric phenomena restricting the availability 
of solar radiation at the earth's surface. The cloudiness as a limiting 
factor causes distribution and dissipation of the solar radiation 
reaching the atmosphere and affects the amount of radiation 
received at the earth's surface. The cloud data are detected 
routinely by meteorological satellites, so a number of models have 
been developed to estimate global solar radiation from observations 
of various cloud layer amounts and cloud types [70]. In the 
following section, correlations which use only the cloudiness are 
presented and classified based on their developing year. 


2.2.1. Model 1: Black model 
Black developed following quadratic equation, using data from 
many parts of the world [71]: 


H 
Ho 


where C is mean total cloud cover during daytime observations 
in octa. 


—0.803—0.340C-0.458C? (C<0.8) (130) 


2.2.2. Model 2: Paltridge model 

Paltridge and Proctor suggested a new model that is able to 
determine the instant and total daily radiation at any location of 
interest. This model takes the solar zenith angle (0), day length 
(So) and cloud factor (CF) as inputs. Their model assumes that the 
effect of atmospheric water vapor, regional albedo and aerosol 
optical air mass on surface radiation is small (less than 5%) [72]. 

The hourly beam (Ip) and diffuse (I4) radiation in (MJ/m? h) are 
determined by: 


Ip = 3.42286[1—exp(—0.075(90—0))] (131a) 


Ia = 0.00913 +0.0125(90—0)+0.723CF (131b) 


The mean monthly Total daily global radiation on horizontal 
surface (MJ/m? day) were computed from correlations below: 


Hi zd EH (1310) 


sunset 


La(Oydt 


sunset 


H=(1-CF) 1,(0) cosOdt + f (131d) 


sunrise sunrise 


In this work, the integral was simply converted to summation 
and added up every 15 min to obtain total daily global radiation. 

The cloud factor (CF) varies from zero for clear sky to 1 for 
overcast sky. This parameter could be obtained by use of the 
numbers of cloudy days in each month and the cloud cover. Cloud 
cover is observed every three hours in three different ranges: (0-2) 
octas, (3-6) octas, and (7-8) octas. To convert the cloud cover to 


cloud factor (CF), the following relationship is used [73]: 


(n1)+4.5(n2)+7.5(n3) 


CF = 
8(n r5 4 n3) 


(131e) 


where nj, n; and n; are the total number of days in each month, 
with zero to 2/8, 3/8 to 6/8, and 7/8 to 8/8 octas, respectively. 


2.2.3. Model 3: Daneshyar model 

Following Paltridge and Proctor’s [72] work, Daneshyar pro- 
posed his method by defining new coefficients for diffuse radia- 
tion adjusted for the climate conditions of Tehran, Iran [74]: 

The hourly beam and diffuse radiation is calculated from 
following equations: 


I, =3.42286[1—exp(-0.075(90-0))] (132a) 


Ia = 0.00515 +.0.00758(90—0) + 0.43677CF (132b) 


The mean monthly total radiation in (MJ/m? day) estimated, 
using time steps of 15 min, is: 


sunset sunset 


H=(1-CP) 1,(0)cosOdt + f Ia(O)dt (1320) 
sunrise sunrise 
2.2.4. Model 4: Badescu model 
Badescu proposed the following correlations [75]: 

H 

Ho a-- bC (133a) 
— =a+bC+cC (133b) 
Ho 

— =a+bC+cC?+dC (1330) 


Ho 


2.2.5. Model 5: Modified Daneshyar model (Sabziparvar model) 
Sabziparvar [76] in 2007 made the following modifications to 
the Daneshyar method [74]: 


(1) In Daneshyar method, solar constant of 1353 (W/m?) has 
been used by the workers. Since the new suggested average 
value of solar constant is about 1367 (W/m?) [77], the total 
daily global radiation was multiplied by a factor. 

(2) For each month, the monthly mean global radiation was 
modified by the Sun-Earth distance correction factor. 

(3) Height effect was also applied separately on beam and diffuse 
radiation. For calculation of the height correction factor in the 
new method, Tehran (where Daneshyar calibrated his sug- 
gested model) is taken as the reference. Details of this method 
are described in [76]. 


2.2.6. Model 6: Modified Paltridge model (Sabziparvar model) 
Similar to method 5, Sabziparvar [76] in 2007 made the same 
modifications to the Paltridge-Proctor method [72]. For locations 
where sunshine duration are not observed, prediction by mod- 
ified Paltridge model which only requires cloud data (are avail- 
able easily by satellites and ground-based measurements) can be 
a good alternative. Details of this method are described in [76]. 


2.3. Temperature-based models 


Cloud observations and sunshine data are not readily available 
in all of the locations. Therefore, developing some precise solar 
radiation models which use commonly available measured para- 
meters such as air temperature is necessary. Due to the common 
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availability of daily maximum and minimum air temperatures, 
several empirical methods have been proposed to estimate solar 
radiation from these variables, especially for locations where the 
air daily temperature is the only available meteorological data. 
The temperature based models assume that the difference in 
maximum and minimum temperature is directly related to the 
fraction of extraterrestrial radiation received at the ground level. 
However, there are factors other than solar radiation that can 
influence the temperature difference such as cloudiness, humid- 
ity, latitude, elevation, topography, or proximity to a large body of 
water [78]. In the following section, the radiation models which 
use the maximum and minimum air temperatures are presented 
and classified based on their years of appearance. 


2.3.1. Model 1: Hargreaves model 

Hargreaves and Samani recommended a simple equation to 
estimate solar radiation using only maximum and minimum 
temperatures [79]: 
H 


g ~” a(Tmax— Tmin)?” 


i (134) 


Initially, coefficient a was set to 0.17 for arid and semi-arid 
regions. Hargreaves [80] later recommended using a=0.16 for 
interior regions and a=0.19 for coastal regions. One of the 
implications of this model for estimating daily global solar 
radiation is: 


- Bayat and Mirlatifi model for Shiraz, Iran [81]: 


H 
I =0.16(Tmax—T min)?” 


a (135) 


2.3.2. Model 2: Bristow and Campbell model 

Bristow and Campbell developed a simple model for daily 
global solar radiation with a different structure in which H is an 
exponential function of AT [82]: 


H 
ms a[1—exp(—bAT*)] (136) 
where AT is the temperature term difference. Although 


coefficients a, b and c are empirical, they have some physical 
meaning. Coefficient a represents the maximum radiation that 
can be expected on a clear day. Coefficients b and c control the 
rate at which a is approached as the temperature difference 
increases. 


2.3.3. Model 3: Allen model 

Following the work of Hargreaves and Samani [79] Allen in 
1997 suggested employing a self-calibrating model to estimate 
mean monthly global solar radiation [78]: 
H 


H = Ki(Tmax—T min)? 


B (137) 


Previously, Allen [83] had expressed the empirical coefficient 
(K,) as a function of the ratio of atmospheric pressure at the site 
(Ps, kPa) and at sea level (Po, 101.3 kPa) as follows: 


P 0.5 
ky = Kra (2) 


For the empirical coefficient k;;, Allen suggested values of 0.17 
for interior regions and 0.20 for coastal regions. Allen reported 
that Eq. (138) performs poorly for sites having an elevation of 
more than 1500 m. 


(138) 


2.3.4. Model 4: Donatelli and Campbell model 

Donatelli and Campbell algorithm is similar in structure to the 
Bristow and Campbell model [82], except taking into account two 
corrective functions of mean and minimum air temperature. The 
model has the following forms [84]: 


H 2 
@ gp a 1cexp(—f (Tave) ATAT min)| (139a) 
Tmax Tmin 
Tavg = Uma 5 ) (139b) 
f (Tavg) = 0.017exp(exp(—0.053T4,4)) (139c) 
f(Tmin) = exp (Tmin/C) (139d) 
H AT" 
(b) m i exp( ba) | (140) 
where AT; is monthly mean AT(°C) 
H C 
(c) Hg a[1—exp(—bf (Tag) AT‘) (141a) 
f (Tag) = 0.017exp(exp(—0.053TagAT*)) (141b) 


2.3.5. Model 5: Hunt et al. model 
Hunt et al. [85] proposed following model by adding another 
coefficient (b) to Hargreaves and Samani model [79]. 


H- Q(Tmax— T min)? Ho +b (142) 


2.3.6. Model 6: Goodin et al. model 

Goodin et al. refined the Bristow and Campbell model [82] by 
adding an Ho-term meant to act as a scaling factor allowing AT to 
accommodate a greater range of H values [86]: 


ac prm) 


The results proved that this model provides reasonably accu- 
rate estimates of irradiance at non-instrumented sites and that 
the model can successfully be used at sites away from the 
calibration site [87]. 


(143) 


2.3.7. Model 7: Thorton and Running model 

The method proposed by Thorton and Running is based on the 
Bristow and Campbell [82] study and is as follows [88]: 
H 


Ho = Tt, maxTf, max (144) 
where Tt, max is the maximum (cloud-free) daily total transmit- 
tance at a location with a given elevation and depends on the 
near-surface water-vapor pressure on a given day of the year, and 
Tf, max Stands for the proportion of Tt, max observed on a given day 
(cloud correction). 


2.3.8. Model 8: Meza and Varas model 

Meza and Varas assumed that a and c coefficients of Bristow- 
Campbell [82] model are fixed and the only b coefficient was 
adjusted to minimize the square errors [89]: 


H 
Ho 


=0.75 [1 exp( bAT’)| (145) 
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2.3.9. Model 9: Winslow model 
The method proposed by Winslow et al. was designated as a 
globally applicable model and the prediction equation is [90]: 


H M 
=1Dl|1-a 
Ho Š | €s(T max) 


(146) 


where es(T min) and es(Tmax) are saturation vapor pressures at min 
and max temperature, respectively. Variable t accounts for atmo- 
spheric transmittance and is estimated from site latitude, eleva- 
tion and mean annual temperature. Function DI corrects the effect 
of site differences in day length, which causes a variation between 
the time of maximum temperature (and minimum humidity) and 
sunset [91]. 
The day-length correction (DI) is approximated by: 


a) í 


day 


Dl = 0 Haay (147) 


2.3.10. Model 10: Weiss et al. model 
Weiss et al. simplified the Donatelli and Campbell [84] and 
Goodin [86] models which only needs one parameter (b) to be 


=a In(Tmax—Tmin)+ b 


ho (153b) 


2.3.14. Model 14: Abraha and Savage model 
Abraha and Savage fixed a=0.75 and c=2 in Donatelli and 
Campbell [84] model and proposed the following correlation [96]: 


H AT? 
Hs =0.75 I zi Oar || 


2.3.15. Model 15: Abraha and Savage, Weiss et al. model 

Weiss et al. and Abraha and Savage put a=0.75 in Donatelli 
and Campbell model [84] and reported the following relationship 
[92,96]: 


i =0.75 [1 exp( bf (Tas) AT^fCT min) ) | 


(154) 


(155) 


2.3.16. Model 16: Almorox et al. model 
A new method is proposed by Almorox et al. for estimating 
daily global solar irradiance in the form of [87]: 


H 


calibrated, while other coefficients (a and c) were fixed as a=0.75 
and c=2 [92]: 


H 

(2g. —075 L exp( bf (Toig) AT?) | (148a) 

f(Tavg) = 0.017exp (exp(—0.053TagAT”) ) (148b) 
H AT? 

(b) Ho =0.75 I zi b Ho )| (149) 


2.3.11. Model 11: Annandale model 

Annandale et al. modified Hargreaves and Samani [79] model by 
introducing a correction factor for parameter a to account the 
effects of reduced altitude and atmospheric thickness on H as [93]: 
A a(1 42.7 x 102) T max—T min)?" (150) 
Ho 
where Z is the elevation above sea level in m. 

Bayat and Mirlatifi used this model with a=0.15 to estimate 
the daily global solar radiation in Shiraz, Iran [81]: 


La 0.15(1 127% 107%) (uncle 


» (151) 


2.3.12. Model 12: Mahmood and Hubbard model 

Mahmood and Hubbard estimated daily solar radiation based 
on maximum and minimum daily air temperatures and proposed 
the following model [94]: 


H- OT max — Tin). Ho” (152a) 
H—2.4999 
Hmod = 0.8023 (152b) 


where Hmoa is estimated solar radiation corrected for systematic 
bias, in MJ/m? day. 


2.3.13. Model 13: Chen et al. model 
Chen et al. presented the following models [95]: 


H 


Ho = A(T max Tmin)” +b 


(153a) 


= A(T max Tand [1 exp( C(€s(Tmin)/@s(Tmax))*) | (156) 


Ho 
where es(Tmin) and es(Tmax) are saturation vapor pressures at 
minimum and maximum temperature, respectively. Saturation 
vapor pressure is related to air temperature and the relationship 
is expressed by [97]: 


e,(T) = 0.61 08exp[17.27(T)/(T +273.3)] (157) 


2.4. Other meteorological parameters-based models 


Accurate prediction of the actual value of solar radiation for a 
given location requires long-term average meteorological data 
which are still scarce especially for underdeveloped and devel- 
oping countries. Therefore it is not always possible to predict the 
solar irradiance for a particular location of interest. Many 
researchers have tried to use various available meteorological 
parameters such as precipitation, relative humidity, dew point 
temperature, soil temperature, evaporation and pressure along- 
side with the classical estimators such as sunshine, air tempera- 
ture and cloudiness to predict the amount of global solar 
radiation. In the following section, the correlations which use 
the various meteorological variables are reported and classified 
based on their developing year. 


2.4.1. Model 1: Swartman and Ogunlade model 

Swartman and Ogunlade stated that the global radiation can 
be expressed as a function of the (S/So) ratio and mean relative 
humidity (RH) [98]: 


HV 

H=a(3,) RH‘ (158a) 
So 

H S 

Ho =a+b($)+cRH (158b) 


2.4.2. Model 2: Sabbagh model 

Sabbagh proposed following equation to estimate the monthly 
average daily global radiation that might be applicable to dry arid 
and semi-arid regions [99]: 


S RH9333 1 


H — 0.06407 (Ks )exp E PE] (159) 
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where Kg is geographical and seasonal factor and can be calcu- 
lated from the following equation: 


0.2 
Ke = 100( 5255750 WpiC03L) 


where Yn is a seasonal factor in month that is suggested by Reddy 
in [100]. 

The model is only reliable for locations with an average mean 
sea level of about 300 m and needs to be modified for arid regions 
with higher altitudes [99]. 


(160) 


2.4.3. Model 3: Gariepy's model 

Gariepy has reported that the empirical coefficients a and b in 
the Angstróm-Prescott [4] model are dependent on mean air 
temperature (T, °C) and the amount of mean precipitation (P, cm) 
[101]: 


a = 0.3791—0.0041 T—0.0176 P (161a) 


b=0.4810+ 0.0043 T+ 0.0097 P (161b) 


2.4.4. Model 4: Garg and Garg model 
Garg and Garg proposed a double linear relation for obtaining 
monthly mean daily global solar radiation as follows [102]: 


H =a+b(=-) +cw 


He s (162a) 


where w(cm) is the atmospheric precipitable water vapor per unit 
volume of air and is computed according to Leckner [103]: 


w= 0.0049RH ES a 3416/ 10] (162b) 
K 
- El-Metwally model for Egypt [7] 
a =0.219+0.526 E +0.004 w (163) 
Ho So 


2.4.5. Model 5: Lewis model 
Lewis reported that global solar radiation on a horizontal 
surface can be calculated by the following equations [104]: 


H=aRH?; logH=a+b logRH (164a) 
H-aS"RH'; logH=a-+b logS+c logRH (164b) 
H =a exp(bS); InH=a+bS (164c) 
H =a exp(bRH); In H= a+bRH (164d) 
H =a expb(S—RH); In H=a+b(S—-RH) (164e) 


2.4.6. Model 6: Ojosu and Komolafe model 
Ojosu and Komolafe proposed the following equation [105]: 


H S Tmi RH 
=a+b EC si) al ) 
Ho (5) (= RH max 


2.4.7. Model 7: Gopinathan model 
Gopinathan introduced a multiple linear regression equation 
of the form [49]: 


(165) 


m - a bcosp- cz «d (S. +eT+fRH (166) 
0 


Ho 


2.4.8. Model 8: De Jong and Stewart model 
De Jong and Stewart introduced the effect of precipitation in a 
multiplicative form as follow [106]: 


H 
Hg = (Tmax Tig) (1 +cP+dP?) (167) 


Hy 


where P is precipitation in mm. 


2.4.9. Model 9: Abdalla model 
Abdalla modified the Gopinathan [49] model for Bahrain as [107]: 


H S 
Ho =a4 (5) E CT +dRH (168a) 
H =a+b + cT -- dRH + ePS (168b) 
Ho So 


Maghrabi [108] estimated the ability of this model for esti- 
mating monthly mean global solar radiation in kW h/m? for 
Tabouk, Saudi Arabia: 


H S 


= —0.107 4- 0.70 0.0025T + 0.004RH (169) 
Ho So 
2.4.10. Model 10: Ododo et al. model 

Ododo et al. proposed two new models as follow [109]: 
H a 
Ho ~ a (=) Tmax RH (170a) 
H S ; S 
uo =e+f| =) +8Tmax+hRH+ iTmax | — (170b) 
Ho So So 


2.4.11. Model 11: Hunt et al. model 

Hunt et al. introduced precipitation (P, mm) and (Tmax, °C) in 
an additive form that have had the highest accuracy at eight sites 
in Canada [85]: 


H= a(Tmax—Tmin) Ho +bTmax+cP+dP? +e (171) 


2.4.12. Model 12: Supit and Van Kappel model 

Supit and Van Kappel proposed a simple method to estimate 
daily global radiation and tested it for various locations in Europe, 
ranging from Finland to Italy [70]. 


H = Hoja y (Tmax—Tmin) + b (1-C/8)]+c 


(172) 


- Supit and Van Kappel model for London weather station, U.K. [70] 


H = Ho[0.061 /(Tmax—Tmin)+0.477,/(1-C/8)]-0.557 (173) 


The regression coefficients of proposed method for various 
locations in Europe are presented in [70]. 


2.4.13. Model 13: Togrul and Onat model 

Togrul and Onat investigated the effect of geographical, 
meteorological and astronomical parameters on the monthly 
mean global solar radiation and gave the following correlations 
[110]: 


H=a+b sinó (174a) 

H=a+bHo (174b) 
S 

H=a+bHo+e(=) (174c) 
0 
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= a+b(5,) +c sinó (174d) 
So 
S = 
H=a+b( 5) +e sind (174e) 
0 
S 
H=a+bHo+c(5,) --dST (174f) 
0 
S 
H=a+bHo4 «(s ) «anu + eST (174g) 
0 
n S 
H=a+b sino+c(5, ) editt cer (174h) 
0 
H=a+bHo+c(5,) -d sind + eRH+fT (174i) 
0 
H=a+bHo+c(5, ) +dRH+esT+fT 174j) 
0 
H=a+bHo+c(5,) +d sind+eRH+fST+gT (174k) 
0 
- Togrul and Onat model for Elazig, Turkey [110] 
H = 4.0891 +6.459 sinó (175a) 
H=-1.32+0.6757Hp (175b) 
H= -1.3876+0.518Hp+2.3064¿) (1750) 
0 
H=2.765+2.2984 (=) +4.9597 sind (175d) 
0 
H=2.6484+4 2.034 (=) + 5.2184 sind—0.0246T (175e) 
0 
H = —1.489 +0.528Ho + 2.537 8 —0.0067ST (175f) 
0 
H = —0.1764+ 0.523Ho + 2.031 (=) —0.0146RH—0.0162ST 
0 
(1758) 
H=3.895+5.116sin ó4- 2.44 (=) —0.0143RH—0.032T 
0 
(175h) 
H= 5.25-0.169H0+2.5017( $) +6.525 sin6—0.0144RH—0.034T 
0 
(175i) 
H = —0.6324-0.516Ho + 2.3255 (=) —0.0103RH 
0 
+0.0373ST—0.062T (175j) 
H = 4.591—0.1135Hpo 4- 2.522 (s. 7- 6.1589 sinó 
0 
—0.0124RH --0.0187ST—0.052T (175k) 


2.4.14. Model 14: Ertekin and Yaldiz model 

Ertekin and Yaldiz estimated the monthly average daily global 
radiation by some multiple linear regression models using nine 
different variables: extra terrestrial radiation, solar declination, 
mean relative humidity, ratio of sunshine duration, mean tempera- 
ture, mean soil temperature, mean cloudiness, mean precipitation 


(cm) and mean evaporation (cm) [2]: 


815 


H=a+bó (176a) 
H=a+bHo+cT (176b) 
S 
H=a+bHo+c(5,) ac (176c) 
0 
S 
H=a+bHo+cRH+a( È) ec (176d) 
0 
S 
H=a+bHo+cRH+d( $) +eT+fC (176e) 
0 
H=a+bHo+cRH+d( $) +eT+fC+iE (176f) 
(0) 
H=a+bHo+có+dRH+e() +fST+gC+hE (176g) 
0 
H=a+bHo+c5-+dRH+e( =) +fT+gST+hC+iE (176h) 
0 
H=a+bHo+c5-+dRH+e( =) +fT+gST+hC+iP+jE (176i) 
0 
- Ertekin and Yaldiz model for Antalya, Turkey [2] 
H=13.58+0.3330 (177a) 
H = —4.46--0.477Hg + 0.226T (177b) 
H = —2.54+ 0.491 Ho + 4.99 (=) —0.406C (1770) 
0 
H = —6.158 +0.487H, + 0.0772RH+ 4.508 (,)-0.6330 (177d) 
0 
H = —7.248 + 0.506Họ + 0.113RH + 5.987 (=) 
0 
—0.075T—0.9204C (177e) 
H = —16.16 + 0.493Ho +0.219RH+9.282 (=) 
0 
—0.247T—0.831C-- 0.2273E (177f) 
H=—11.90+0.353H 9 + 0.1096 +0.229RH 
+11.095 (=) —0.279ST—0.995C+0.241E (1778) 
0 


H = —12.57 +0.361Ho +0.09976+0.234RH+ 11.19 (5 


So 
0.069T—0.227ST—0.99C 4- 0.251E 


H = —13.08 + 0.386Ho + 0.0902 + 0.2254RH 4- 11.59 6 


0 


—0.034T + 0.251ST—0.977C—0.0072P+ 0.2373E 


s) 


) 


(177b) 


(177i) 


Menges et al. calibrated Eq. (176i) for Konya, Turkey and 


obtained the empirical coefficients of this model as follow 


H = 20.296019—0.096134Hp + 0.3175936 


—0.146422RH 4- 10.705159 (5) 
0 


—0.288332T +0.021331ST+0.359791C 
+0.207588P—0.076444E 


[111]: 


(178) 
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2.4.15. Model 15: Trabea and Shaltout model 
Trabea and Shaltout introduced the following correlation to cal- 
culate the daily global solar radiation at five stations in Egypt [112]: 


ff: =a+b 2 +CT max +dV+eRH+fPS (179) 
Ho So 
- Trabea and Shaltout model for Egypt [112] 
H = —0.1394- 0.229 3 +0.009T max + 0.004 V 
Ho So 
+0.002RH + 0.002PS (180) 


2.4.16. Model 16: El-Metwally model 

El-Metwally proposed following simple new methods to esti- 
mate global solar radiation based on meteorological data over six 
stations in Egypt [113]: 


H = aHo+bT max +CT min + dC +e (181a) 
H = aHo + bT max + CT min + dC (181b) 
H = exp(aHo + bT max +CT min + dC +e) (181c) 


2.4.17. Model 17: Chen et al. model 

Chen et al. suggested a logarithmic relationship between the 
daily global solar radiation, daily extraterrestrial solar radiation 
and the temperature difference between the maximum and 
minimum daily air temperature as follow [95]: 
H 


sy 
— =a In(Tmax—T min) +b (=) +d 
0 


a (182) 


2.4.18. Model 18: Modified Sabbagh method (Sabziparvar model) 
Sabziparvar in 2007 revised Sabbagh [99] model to predict the 
monthly average daily solar radiation on horizontal surfaces in 
various cities in central arid deserts of Iran. The modification was 
made by the inclusion of following factors to Sabbagh model: 
(1) Height correction factor, (2) Sun-Earth distance correction 
factor and (3) inclusion of monthly total number of dusty days [76]. 


2.4.19. Model 19: Bulut and Büyükalaca model 

Bulut and Büyükalaca proposed a simple model for estimating 
the daily global radiation on a horizontal surface. The model is 
based on a trigonometric function, which has only one indepen- 
dent parameter, namely the day of the year [114]: 


T 
365 


They examined their model for 68 provinces of Turkey [114]. 
Some investigator employed this model to estimate global 
solar radiation in different locations as follows: 


I=b+(a b)|sin| ees] (183) 


- Bulut and Büyükalaca model for Ankara, Turkey [114] 


T 


1.5 
I= 3.86+(22.71-3.86)|sin sgg "+ 5)| | (184) 
— El-Sebaii et al. model For Jeddah, Saudi Arabia [19] 
4.92+9.61|sin|.~ 5 8 
1—14924-9. sin [362 (1+5)]| (185) 


2.4.20. Model 20: El-Sebaii et al. model 
El-Sebaii et al. developed and presented empirical correlations 
between the monthly average of daily global solar radiation 


fraction (H/Ho) and various meteorological parameters [19]: 


H S 
Fy et (3) «er (186a) 
H S 
Ho = a+b(5,) +cRH (186b) 
H 
+ =a+bT+cRH (186c) 
Ho 
H 
Ho = a+ b(Tmax—Tmin) + CC (186d) 
H , 05, 

=a+D(Tmax—Tmin) ^ 4 CC (186e) 
Ho 
H S 
Fe et (s) «ec (186f) 


— El-Sebaii et al. model For Jeddah, Saudi Arabia [19] 


H 1924260 (=) +0.006T (187a) 

Ho So 

H s 

= 1.624224( >) +0.332RH (187b) 

Ho So 

a =0.139-0.003T+0.896RH (1870) 
0 

H 

gj; = 0.214+0.035(Tmax—Tmin)—0.028C (187d) 
0 

H 2405 

E = 0.08 + 0.21(Tmax—Tnin)°®—0.012C (1876) 
0 

H .-2764372( 2) 40.0010 (187f) 

Ho So 


2.4.21. Model 21: Maghrabi model 
Maghrabi established a simple model to calculate the monthly 
mean global solar radiation on a horizontal surface using five 
meteorological parameters [108]: 
H=a+b($ )+cT+dPo+e PWV+fRH (188) 
0 


The precipitable water vapor (PWV) is calculated using the 
Reitan [115] equation as follows: 


In PWV (mm) = 0.1102 +0.0613 ty (189) 
- Maghrabi model for Tabouk, Saudi Arabia [108] 
H(kw h/m?) = 163.01 1.04( =") +0.12T—0.21Po 
0 
—1.06PWV —0.03RH (190) 


3. A case study 
3.1. Experimental data 


It is already mentioned that the main objective of this study is 
to comprehensively collect and review the global solar radiation 
models available in the literature and categorize them based on 
the employed meteorological parameters. In order to evaluate the 
applicability and accuracy of the collected models for computing 
the monthly average daily global solar radiation on a horizontal 
surface, the geographical and meteorological data of Yazd, Iran 
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was used. Yazd city is located in latitude 31°54’N, longitude 
54^17'E and its altitude of sea level is 1237.2 m. Yazd province 
is situated in the central part of Iranian plateau and because of 
being in the sun-belt, it is an ideal location to benefit the 
advantages of solar energy utilization and adoption of its related 
technologies [116]. The total radiation received by a horizontal 
surface in Yazd City during a year is reported to be around 
7787 Mj[m? with about 3270 sunshine hours in a year while the 
number of overcast days is fewer than 1110 h per year [117]. In 
the present work, the value of measured data of daily global solar 
radiation on a horizontal surface in Yazd city in the period of 
1982-2008, as well as number of sunshine hours, maximum, 
minimum and daily ambient temperatures, amount of cloud 
cover, relative humidity and the other related meteorological 
parameters during the period of 1988-2008 were obtained from 
the Islamic Republic of Iran Meteorological Office (IRIMO) data 
centre. 

The measured global solar radiation data were checked and 
controlled for errors and inconsistencies. The purpose of data 
quality control was to eliminate faulty data and inaccurate 
measurements. After the quality control, the measured data 
were averaged to obtain the monthly mean daily values by taking 
the data for the average day of each month as recommended 
by Duffie and Beckman [118]. The measured data were then 
divided into two sets. The first sub-data set (1988-2003) were 
employed to develop empirical correlations between the monthly 
average daily global solar radiation fraction (H/Ho) and monthly 
average of desired meteorological parameters (for the global solar 
radiation this subset is contained the measured data of 1982- 
2003). The second sub-data set (2004-2008) were then used to 
validate and evaluate the derived models and correlations. 
The measured long-term monthly average daily global radiation 
distribution for Yazd city in the periods of 1982-2008 is shown 
in Fig. 1. 


3.2. Statistical evaluation 


The accuracy and performance of the derived correlations in 
predicting of global solar radiation was evaluated on the basis of 
the following statistical error tests which are coefficient of 
determination (R?), root mean square error (RMSE), mean bias 
error (MBE), mean absolute bias error (MABE), mean percentage 
error (MPE), correlation coefficient (r) and t-Test statistic. These 
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Fig. 1. The measured values of monthly average daily global solar radiation for 
Yazd city in the periods of 1982-2008. 


error indices are defined as: 


Yi- 1 (Him —Hic)” 


R=1 5 (191) 
Viet (Him-Hm) 
1 1/2 
RMSE= By (Ha)! (192) 
i= 
1 n 
MBE = D» (Hic—Him) (193) 
1 n 
MABE= 124 | (Hic—Him) | (194) 
124 /H;,—H; 
MPE(%) = — Mein) 100 195 
w=) ( n) x (195) 
iat (Hi,.—Hc) (Him —Hm) (196) 


Je b» m (Hi, —Hc)* Xii (Him—Hm) 14 


In the above relations, the subscript i refers to the ith value of 
the solar irradiation and n is the number of the solar irradiation 
data. The subscripts “c” and “m” refer to the calculated and 
measured global solar irradiation values, respectively. H¿ is the 
mean calculated global radiation and Hm is the mean measured 
global radiation. 

Although these statistical indicators generally provide reason- 
able criteria to compare models, but do not objectively indicate 
whether the estimates from a model are statistically significant. 
The t-statistic (Eq. (197)) allows models to be compared and at 
the same time indicates whether a model's estimate is statisti- 
cally significant at a particular confidence level or not. 


2 112 
-| (n—1)(MBE) | a97 


(RMSE)? — (MBE)? 


The smaller the t-value, the better the model performance is. 
To determine whether the estimates from a model are statistically 
significant, one simply has to determine from the standard 
statistical tables, the critical t-value, i.e., (t,/2) at « level of 
significance and (n—1)° of freedom. The model is judged to be 
statistically significant if the calculated t-value is less than the 
critical value. 


3.3. Result and discussion 


3.3.1. Calibration 

Linear, multiple linear and nonlinear regressions were carried 
out between the monthly mean measured global solar radiation 
and the meteorological parameters using the first sub-data set to 
obtain the values of empirical coefficients of the selected models 
from each category. The coefficient of determination (R?) index is 
used to determine that how well the regression line approximates 
the real data points. A model is more efficient when R? is closer to 
1. In addition, the t-value index is used to determine whether the 
estimates from a model are statistically significant or not. For the 
model's estimates to be judged statistically significant at the 
(1—«) confidence level, the calculated t-value must be less than 
the critical value. Regression coefficients of the selected models in 
each category for Yazd city along with the values of coefficients of 
determination (R?) and t-values are presented in Table 1. It should 
be mentioned that several radiation models from each category is 
used for estimating the monthly average daily global solar 
radiation fraction (H/Ho), however, the results of one selected 
model from each category is presented in Table 1 for brevity. 
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Table 1 


Regression coefficients of selected radiation models for Yazd city along with the values of coefficients of determination (R?) and t-value. 


Models (R2) t-value a b c d 
Sunshine-based models 
El-Metwally model (Eq. (116)) 0.979 0.5234 0.7272 
Cloud-based models 
Badescu model (Eq. (133a)) 0.9948 0.2663 0.7658 — 0.0594 
Temperature-based models 
Hargreaves model (Eq. (134)) 0.9543 0.2625 0.1746 
Other meteorological parameters-based models 
Chen et al. model (Eq. (182)) 0.9971 0.3782 — 0.01088 5.127 0.101 — 4.299 


Table 2 
Statistical results for the validation of the selected models for Yazd city (using data 
in the period of 2004-2008). 


Models RMSE MBE MABE MPE r 
(MJ/m? (MJ/m? (MJ/m? (%) 
day) day) day) 


Sunshine-based models 


El-Metwally model 0.5385 0.04371 0.4124 0.8485 0.9969 
Cloud-based models 

Badescu model (a) 1.152 —0.3385 0.9592  —1.877 0.9846 
Temperature-based models 

Hargreaves model 0.7103 —0.02285 0.6293 0.3039 0.9934 
Other meteorological parameters-based models 

Chen et al. model 0.8542 0.2104 0.746 0.9859 0.9913 

35 
—— H measured 


mM El-Metwally model [7] 
Xx Badescu model [75] 
A Hargreaves model [79] 


e Chen et al. model [95] 


H (MJ/m2day) 


Jun. 4 
Jul. 

Aug. 4 
Nov. + 


May. 4 


a 
< 


Jan 
Feb. 4 
Mar. 4 


g 

o 

a 
Month 


Fig. 2. Measured (in the periods of 2004-2008) and calculated values of monthly 
average daily global solar radiation for Yazd city. 


The results of Table 1 show that all the selected models give 
the coefficients of determination (R?) higher than 95%, indicating 
a very good fitting between the monthly average daily global 
radiation and the other meteorological parameters. The best fit is 
obtained by employing Chen et al. model [95]. Comparison 
between the different models according to the t-value shows that 
the calculated t-values are less than the critical t-value (2.201 at 
95% confidence level), indicating that all models have statistical 
significance. 


3.3.2. Validation 

The performance of developed correlations in predicting the 
global solar radiation on horizontal surface is also assessed by 
comparing each models outputs and the measured values from 
the second subset data based on the r, RMSE, MBE, MABE and MPE 
error indices. A summary of all the statistical parameters are 
presented in Table 2. For higher modeling accuracy RMSE, MBE, 


MABE and MPE indices should be closer to zero, but correlation 
coefficient (r) should approach to 1 as closely as possible. 

The performance of the most accurate models in each category 
for estimating the monthly average daily global solar radiation in 
Yazd is presented and compared in Table 2. From the statistical 
analysis, it can be seen that the estimated values of monthly 
mean daily global solar irradiation are in good agreement with 
the measured values for all models. It is found that the values of 
RMSE are in the range from 0.5385 to 1.152 (in MJ/m? day). The 
MBE achieved in this study for all models are in the acceptable 
range. Negative values of MBE in temperature and cloud based 
models indicate an underestimation of measured global solar 
radiation by these models. As an over-estimation of an individual 
observation may cancel under-estimation in a separate observa- 
tion, using MABE index is more appropriate than MBE. The biggest 
value of MABE is belonged to Badescu model (a) with 0.9592. The 
mean percentage errors (MPE) of all models are in the range of 
acceptable values between 0.3039% and — 1.877%. Also according 
to the statistical test of the correlation coefficient (r), all models 
give very good results (above 0.98). El-Metwally sunshine based 
model [7] gives the highest value of r with value of 0.9969, 
indicating a strong positive linear relationship between the 
measured and calculated values of the global solar radiation. 

The measured values of the monthly average daily global 
radiation and corresponding calculated values by employing 
El-Metwally model [7], Badescu model [75], Hargreaves model 
[79] as well as Chen et al. model [95] for Yazd city are illustrated 
in Fig. 2. As may be seen, agreement between the values obtained 
from the selected models and the measured data are very good. 

Based on the statistical indicators presented in Table 1 as well 
as Fig. 2, all models show good estimation of the monthly average 
daily global solar radiation on a horizontal surface for Yazd city. 
But the El-Metwally sunshine based model [7] produces the best 
result among the correlations developed in this study. Therefore, 
this model is recommended for the prediction of global solar 
radiation on a horizontal surface in Yazd and elsewhere with 
similar climatic conditions, where the radiation data are missing 
or unavailable. 

The agreement between measured and calculated value of H 
was also confirmed for each month by calculating the relative 
percentage error (e). The relative percentage error for each month 
is defined as: 


(198) 


Fig. 3 shows the comparisons of the relative percentage error 
(e) of measured and calculated H for all developed models. It is 
seen from Fig. 3 that the relative percentage error for each month 
rarely exceeds + 10%. 

It is found from Fig. 3 that almost all methods provide low 
performance at both winter and spring, this is due to the increase 
of cloudiness and aerosol content, respectively. 
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Fig. 3. Comparisons of monthly relative percentage error (e) of monthly average 
daily measured and calculated global solar radiation for Yazd (2004-2008). 


4. Concluding remarks 


Solar energy technologies offer opportunities for use of a clean, 
renewable, and domestic energy resource, and are essential 
components of a sustainable energy future. Availability of a 
complete and accurate solar radiation data base for each specific 
region is essential in the prediction, study, and design of solar 
energy systems. Information on global solar radiation at any sites 
(preferably obtained over a long period of time) is useful not only 
to the locality where the radiation data is collected but also for 
the wider world community. 

A global study of the world distribution of global solar 
radiation requires knowledge of the radiation data in various 
countries. For the purpose of worldwide marketing, the designers 
and manufactures of solar equipment will need to know the 
average global solar radiation available in different regions. In this 
regard, a common practice is to estimate average daily global 
solar radiation using appropriate correlations which are empiri- 
cally established by using the measured data at some selected 
locations. These correlations estimate the values of global solar 
radiation from more readily available meteorological, climatolo- 
gical, and geographical parameters. The number of such correla- 
tions or empirical equations that have been developed and tested 
is relatively high, which makes it difficult to choose the most 
appropriate correlation for a particular purpose and site. 

In this study, a comprehensive collection and review of available 
solar models is conducted to assist in the selection of most 
appropriate and accurate model based on the available measured 
meteorological data. Moreover, the collected models are classified 
into four categories namely, sunshine-based, cloud-based, tempera- 
ture-based, and other meteorological parameter-based models. The 
regression constants of some collected solar models have been 
generally presented to calculate the global solar radiation with high 
accuracy in a given location. These regression coefficients have been 
computed using measured monthly mean daily global radiation and 
other meteorological variables for locations of interest. The regres- 
sion coefficients in the presented models differ from location to 
location. This may be due to the local and seasonal changes in the 
type and thickness of cloud cover, the effects of snow covered 
surfaces, the concentrations of pollutants, and latitude [13]. 

The correlations reported in this study have had the best 
estimation in the location of interest and the present classification 
of solar models make it easy to choose a suitable model from 
various existing models based on the meteorological parameters 
measured in the location of interest. These correlations with these 
selected variables can be used by solar energy system designers in 
the design and performance prediction of solar applications with 
fairly high level of accuracy in calibration sites and possibly in 
regions with similar climatic conditions. 


Furthermore, to evaluate the accuracy and applicability of 
collected various models for estimating the monthly average 
daily global radiation on a horizontal surface, the long term 
measured meteorological and solar data of Yazd city, Iran is 
employed. With employing linear and nonlinear least square 
method, regression coefficients of some selected models from 
the four categories were obtained by using first sub-data set. The 
developed models were then compared with each other and with 
the experimental data of second subset on the basis of the 
statistical error indicators such as RMSE, MBE, MABE, MPE and 
correlation coefficient (r). The most accurate models in each 
category are then identified and presented. These models have 
reasonable values of estimation errors. Based on the statistical 
results, El-Metwally sunshine-based model reproduce the mea- 
sured monthly average daily global solar irradiation for Yazd city 
with the highest accuracy. 

A common feature of almost all models collected in this paper 
is that they account for latitude, solar declination, elevation, day 
length and atmospheric transmissivity by including the extra- 
terrestrial radiation (Ho) term in the model. Among the four 
categories of solar radiation models, the sunshine-based methods 
are generally more accurate but they are often limited by the lack 
of availability of sunshine records. Estimation of global solar 
radiation from the air temperature offers an important alternative 
in the absence of measured H or sunshine duration because of the 
wide availability of air temperature data. The main advantage of 
models in this category is the readily available data. Temperature- 
based models are a convenient tool for estimating solar radiation 
if the parameters can be calibrated for each specific location [87]. 
The temperature range (AT) is the main factor affecting accuracy 
of the temperature-based models. Larger (AT) generally results in 
a better predictive accuracy, meaning that the temperature-based 
models are more applicable in areas with larger temperature 
range. This implies that model calibration is particularly sensitive 
in humid regions where (AT) is generally small [119]. 

When sunshine duration data are not available, prediction by 
the method that only requires cloud data (are available easily by 
satellites and ground-based measurements) can be a good alter- 
native. One drawback of these models is that they are sensible to 
human biasing [70]. 

The models presented in the last category, employ more than 
one meteorological parameter for predicting the global solar 
radiation. These models reported a good estimation of solar 
radiation, but they are limited because of using the various 
meteorological parameters as models inputs that are not readily 
available in most of the location of interest. Moreover, to prevent 
the error of measurement tools also increasing speed, the solar 
radiation should be estimated using the minimum measured 
input parameters. 
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